Subdominant matroid ultrametrics 

Federico Ardila 
Abstract 

Given a matroid M on the ground set E, the Bergman fan B{M)^ or 
space of A/-ultrametrics, is a polyhedral complex in which arises in 
several different areas, such as tropical algebraic geometry, dynamical 
systems, and phylogenetics. Motivated by the phylogenetic situation, 
we study the following problem: Given a point uj in M^, we wish to 
find an M-ultrametric which is closest to it in the ^oo-metric. 

The solution to this problem follows easily from the existence of the 
subdominant M-ultrametric: a componentwise maximum M-ultrametric 
which is componentwise smaller than ld. A procedure for computing 
it is given, which brings together the points of view of matroid theory 
and tropical geometry. 

When the matroid in question is the graphical matroid of the com- 
plete graph Km the Bergman fan B{Kn) parameterizes the equidistant 
phylogenetic trees with n leaves. In this case, our results provide a con- 
ceptual explanation for Chepoi and Fichet's method for computing the 
tree that most closely matches measured data. 

1 Introduction 

Given a matroid M on the ground set E, the Bergman fan B{M), or space 
of M-ultrametrics, is a polyhedral complex in which arises in several 
different areas, such as tropical algebraic geometry jj^l) dynamical systems 
|H] and phylogenetics p. It has been described topologically and combina- 
torially pQ. Motivated by the phylogenetic situation, we study the following 
problem: Given a point lv in M^, we wish to find an M-ultrametric which 
is closest to it in the ^oo-metric. 

In SectionEJwe define the Bergman fan B{M), as well as the notion of an 
M-ultrametric. We offer several characterizations of them. When M{Kn) 
is the graphical matroid of the complete graph Kn, M(i^„)-ultrametrics are 
precisely ultrametrics in the usual sense. 
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In Section El we show that, in fact, there is a componentwise maximum 
M-ultrametric which is componentwise smaller than uj. We call it the suh- 
dominant M-ultrametric of u, and denote it u!^ . A simple translation of it 
is a closest M-ultrametric to to. A procedure for computing subdominant 
M-ultrametrics is given, similar in spirit to Tarjan and Kozen's blue and 
red rules ^ for computing the bases of minimum weight of a matroid. 

In Section 13 we prove that the Bergman fan is a tropical polytope in the 
sense of [7]j and that the subdominant M-ultrametric o;*^ is precisely the 
tropical projection of uj onto B{M). 

In Section |SJ we discuss a special case of particular importance: the 
Bergman fan of M{Kn), the graphical matroid of the complete graph Kn- 
As shown by Ardila and Klivans (Ij, the Bergman fan B{Kn) can be re- 
garded as a space of phylogenetic trees, and we are interested in finding an 
(equidistant) phylogenetic tree that most closely matches measured data. 
In this case, our results provide a conceptual explanation for the method 
developed by Chepoi and Fichet in fB; to compute the tree that most closely 
matches measured data. 

Throughout this paper, familiarity with the fundamental notions of ma- 
troid theory will be very useful; we refer the reader to pTO:, Ch. 1,2] for the 
basic definitions. 

2 The Bergman fan and matroid ultrametrics 

Let M be a matroid of rank r on the ground set E, and let u> € M.^ . Regard 
a; as a weight function on M, so that the weight of a basis B = {bi, . . . ,br} 

of M is given by = + "I \-i^br- 

Let M(j be the collection of bases of M having minimum w-weight. This 
collection is itself the set of bases of a matroid; for more information, see 
for example ^j. 

Definition. The Bergman fan of a matroid M with ground set E is: 

I3{M) := {cj E : every element of E is in a uo-minimum basis}. 
An M-ultrametric is a vector in B{M). 

Our ongoing example throughout the paper will be the Bergman fan of 
the matroid M^K/^), the graphical matroid of the complete graph K/j^. An 
M(/C4)-ultrametric is an assignment of weights to the edges of K4 such that 
any edge of -fC4 is in a spanning tree of minimum weight. An example of an 
M(i('4)-ultrametric is given in Figure^ The w-minimum spanning trees are 



2 



A 




B 0.2 c 



Figure 1: An M(i^4)-ultrametric. 



those consisting of an edge of weight 2, an edge of weight 1.2, and an edge 
of weight 0.2. Every edge of the graph is in at least one such tree. 

We start by reviewing some useful facts about matroids. Given a basis 
B and an element x ^ B, there is a unique circuit which is contained in 
BUx (and must contain x). It is called the fundamental circuit of B and x. 

Given a basis B and an element y (z B, there is a unique cocircuit which 
is disjoint with B — y (and must contain y). It is called the fundamental 
cocircuit of B and y. It is equal to the fundamental circuit E — B and y 
in the dual matroid M* . 

Proposition 1. /5, Lemma 7.3.1] Let B he a basis of M , and let x ^ B 

and y & B. The following are equivalent. 

(i) B\J X — y is a basis. 

(a) y is in the fundamental circuit of B and x. 
(Hi) X is in the fundamental cocircuit of B and y. 

We will also need the following lemma: 

Lemma. 'JO, Prop. 2.1.11] A circuit and a cocircuit cannot intersect in 
exactly one element. 

We now give two additional characterizations of the Bergman fan, which 
will be central to our analysis. 

Proposition 2. Let M be a matroid with ground set E, and let lo € 'Mp . 

The following are equivalent: 

(i) uj is an M -ultrametric. 
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(ii) No circuit has a unique uj-maximum element. 
(Hi) Every element of E is uo-minimum in some cocircuit. 

Proof, {ii) (i): Let e £ E; we want to show that it is in an w-minimal 
basis. Consider any w-minimal basis B. If e G B, there is nothing to prove. 
Otherwise, the fundamental circuit of B and e has at least two w-maximum 
elements; let / 7^ e be one of them. By Proposition^ i3 U e — / is a basis; 
since ujf > uJe, it is actually the cj-minimum basis that we need. 

(Hi) =^ (ii): Assume that a circuit C achieves its w-maximum only at e. 
Say C* is a cocircuit where e is w-minimum. By the Lemma, we can find an 
element / 7^ e in C n C*; then uje> uJf because e is the unique w-maximum 
in C, and uJe <: uJf because e is cj-minimum in C*. 

(i) =^ (Hi): Let e G £"; we want to show that it is cj-minimum in some 
cocircuit. Let B be an w-minimum basis containing e, and let / be the 
cj-minimum element in the fundamental cocircuit C* of B and e. Then 
U / — e is a basis, and its weight is at least ljb- Therefore tof > tOe, so e 
is w-minimum in C* also. 

□ 

Let us check Proposition[51for the M(i^4)-ultrametric of Figure^ State- 
ment (ii) is clear: in each cycle of K4, the two largest weights are equal. To 
check statement (Hi), recall that the cocircuits of M{K4) are the cuts of K4. 
Denote hy S — S' the cut that separates the vertices in S from the vertices in 
S' . Then the edges of weight 2 are minimum in the cut ABC — D, the edges 
of weight 1.2 are minimum in the cut A — BCD, and the edge of weight 0.2 
is minimum in the cut AB — CD. 

3 Subdominant M-ultrametrics. 

For uj, uj' G M^, say that w < if < u)'^ for each e £ E. 

Proposition. Let M be a matroid on E, and let uj G M^. There exists a 
unique maximum M -ultrametric which is less than or equal to to. We call it 
the subdominant M-ultrametric of oj, and denote it a;*^. 

Proof. Let S = {uj' G B{M) : uj' < u>}, and let = sup 5, where the sup 
is taken componentwise. We claim that uj^'^ G B{M). 

Proceed by contradiction; assume that e is the unique cj^^-maximum of 
the circuit C. Let e be such that w*^ — e > uj^ for all / G C — e. We can find 
a u;' G S such that > - e. But then > u:¥ > uj'. for ah / G C - e. 
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so e is the unique w'-maximum of C. This contradicts the assumption that 
uj' G B{M). □ 



We now turn to the problem of constructing the subdominant M-ultra- 
metric of a vector oo. Take a vector uj which is not an M-ultrametric. This 
means that conditions {ii) and [iii) of Proposition [21 are not satisfied. 

We must force each element to be w-minimum in some cocircuit, and 
not allow it to be the unique w-maximum element of a circuit. At the same 
time, we must demand that it keeps its weight as high as possible. Let us 
issue the following two rules, similar in spirit to Tarjan's blue and red rules 
for constructing the minimum weight spanning trees of a graph. [HI 

Blue Rule: Suppose an element e is not w-minimum in any cocircuit. 
Look at the minimum weight in each cocircuit containing e. (They are 
all less than uje-) Make LOe equal to the largest such weight. 

Red Rule: Suppose an element e is the unique w-maximum element 
of a circuit. Look at the maximum weight in C — e for each circuit C 
containing e. (At least one of them is less than coe-) Make equal to 
the smallest such weight. 

Figure |21 shows an example of an assignment uj of weights to the edges 
of which is not an ultrametric, and the result of applying the blue rule 
or the red rule to edge CD. When the blue rule is applied, the edge CD 
inherits its new weight from the cut ABC — D, as shown. When the red rule 
is applied, the edge CD inherits its new weight from either one of the cycles 
ACD (as shown) or ABCD. Notice that, surprisingly, the blue rule and the 
red rule give the same result. The reason for this will soon be explained. 
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Figure 2: Applying the blue rule or the red rule to edge CD. 

In principle, it is not yet clear how to apply these rules, or what they 
will do. We do not know in what order we should apply them, and different 
orders would seem likely to give different results. It is not even obvious that 
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applying these two local rules successively will accomplish the global goal of 
turning uj into an M-ultrametric. 

Fortunately, as the following theorem shows, the situation is much sim- 
pler than expected. It turns out that, if each element of E is willing to do 
its part by complying with the two rules, the global goal of constructing the 
subdominant M-ultrametric will inevitably be achieved. 

Theorem 1. Let M be a matroid on the ground set E, let uo G . For 

each element e (z E, there are two possibilities: 

1. Both the blue rule and the red rule apply to e, and they both change 
its weight from to w*^, or 

2. Neither the blue rule nor the red rule apply to e, and oje = uj^ . 

Consequently, if we apply the blue rule or the red rule to each element of 
E in any order, we obtain the subdominant M-ultrametric lj^P 

Before proving Theorem ^ let us illustrate it with an example: the 
construction of the subdominant ultrametric of the weight vector to we con- 
sidered in Figure l2j Figure 01 shows the result of applying, at each step, 
either the blue rule or the red rule to the highlighted edge. The reader is 
invited to check that, at each step of the process, the blue rule and the red 
rule give the same result. 
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Figure 3: Constructing the subdominant M(i('4)-ultrametric of u;. 

After three steps, we reach an ultrametric. This means that we are done: 
we have in fact reached u!^^^^\ the subdominant M(i<'4)-ultrametric of lv. 
Moreover, we could have applied either the blue rule or the red rule to the 
edges of K4 in any order, and we would have obtained exactly the same 
result: uj^^^^-^l 

Proof of Theorem^ Let lo'^ and w'^ be the weights assigned to e by the blue 
rule and the red rule, respectively. If the blue rule (or the red rule) does not 
apply to e, set = Wg (or cjg = uig). We proceed in several steps. 
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1. uj>uj' > LO^"^ and uj>u)"> uj^ . 

The inequalities > cjg > uj^ and uJe > oj'l > oj^ are easy to see: 
When one of the rules applies, it decrees a decrease in the weight of e that 
is clearly necessary for any M-ultrametric which is less than uo. 

2. u' is an ultrametric. 

We use Proposition l^fiu). Let e & E. By definition, is equal to the 
w-minimum weight of some cocircuit C* containing e. (This is true even 
if the blue rule does not apply.) We claim that it is also the w'-minimum 
weight of C*. Let f (z C*. Then iv'j is the largest w-minimum weight of 
a cocircuit containing /; this is greater than or equal to the w-minimum 
weight of C*, which is Wg. 

3. oj' = iv^^. 

We already knew that u' > lo^ . Now, since u' is an ultrametric less 
than or equal to u>, it must be less than or equal to uj^ by Proposition |31 

4. to" = io^. 

We already knew that u" > uj^'^ . Now assume that > t > oj'^. Let 
S = {f & E\ujf < t}. We have Wg > uj" > t, so e is not in S. Since 
Wg > t, no circuit containing e is in S" U e. Therefore e ^ cl(S'), and there 
is a hyperplane H containing cl{S) with e ^ H. The complement of iJ is a 
cocircuit C* which contains e, and does not contain any element of S; its 
w-minimum weight is at least t. Therefore uj'^> t, a contradiction. □ 

If we have an w-minimum basis of M, which is easy to construct using 
Tarjan and Kozen's blue and red rules |14| l9]. the task of constructing uj'^^ 
becomes even simpler, as the following result shows. 

Proposition. Let B he any tu-minimum basis. 

(i) For each e G B, lo^^ = uj^- 

(a) For each e ^ B, let Cq be the fundamental circuit of B and e. Then 
Wg^ is equal to the uj-maximum weight in Cq — e. 

Proof. It follows from Proposition ^ that any element e G i? is w-minimum 
in the fundamental cocircuit of B and e. Therefore the blue rule does not 
apply to e, and (i) follows from Theorem^ 

Now we show that, for e ^ B, miuggc" maxjgc'_e loj = maxj^Co-e ^j- 
Say the maximum in the right hand side is achieved at uJm- First we claim 
that, for any circuit C containing e, there is a basis B — mUf with / G C — e. 
In fact, B — mU C has full rank, since it contains the basis B — mU e. But 
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e is dependent on C — e, and hence on {B — m) U (C — e) also. Therefore 
{B — m) U (C — e) has full rank as well, and the claim follows. Now, since 
S — m U / is a basis, < loj < maxj^c-e ^j- The desired equality follows. 

Now (ii) follows easily. If the red rule applies to e, it will change iVe 
to LO^^ = oJm. If it does not apply, then e is not the unique w-max of any 
circuit, and co^'^ = uJe = iOm- □ 

4 Tropical projection 

In this section we study the Bergman complex of a matroid from the point 
of view of tropical geometry. We start by reviewing the basic definitions; 
for more information, we refer the reader to |7j. 

The tropical semiring (M U {oo},0,0) is the set of real numbers aug- 
mented by infinity, together with the operations of tropical addition and 
multiplication, which are defined by x © y = min(x, y) and x Q y = x + y. 

Throughout this section, the symbol M will denote the tropical semiring, 
and not the ring of real numbers. The additive and multiplicative units of 
M are oo and 0, respectively. 

The n-dimensional space M" is a semimodule over the tropical semiring, 
with addition x (B y and scalar multiplication cQ x given componentwise. 
Tropical scalar multiplication by c is equivalent to translation by the vector 
(c, c, . . . , c). The (n — 1)- dimensional tropical projective space is TP"^"'^ = 
W/{x ~ (cQx)). 

A subset S of is tropically convex if (a©x)®(60y) € S for any x,y € S 
and any a,b G M. Tropically convex sets are invariant under tropical scalar 
multiplication; i.e., translation by (c, c, ...,c). We will therefore identify 
them with their image in TP"^"*^. A tropical polytope is the set 

tconv(y) = {(ai Qvi)®---e{arQVr) : oi, . . . , o,. G M} 

of all linear combinations of a finite set V = {vi, . . . ,Vr} ^ TP"~^. 

Our claim is that, for any matroid M on the set [n], —B{M) is a tropical 
polytope in TP""^ For each flat F of M, let vp G TP""^ be the vector whose 
z-th coordinate is oo if i G -F, and otherwise. Recall that the hyperplanes 
of a matroid are its maximal proper flats. 

Proposition. For any matroid M on [n], —B{M) is a tropical polytope in 
']pp«-i_ Jig of vertices is 

Vm = {vh : H is a hyperplane of M}. 
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Proof. We start by reviewing the description of the Bergman fan obtained 
in [T]. Given a flag of subsets = {0 =: Fq C Fi C • • • C C F^+i := E}, 
the weight class of is the set of lij G for which uj is constant on each set 
Fi — Fj_i, and has uj\f^-f^_i < uj\F^+i-Fi- For example, one of the weight 
classes in is the set of vectors uj such that loi = L04 < uj2 < ^3 = ^5- It 
corresponds to the flag {0 C {1, 4} C {1, 2, 4} C {1, 2, 3, 4, 5}}. 

The disjoint union of all weight classes is M". The Bergman fan B{M) 
is also a union of disjoint weight classes: the weight class of ^ is in ;B(M) if 
and only if ^ is a flag of flats of M. 

Now we will show that 

—B{M) = tconvl^j? : F is a proper flat of M}. 

The left hand side is contained in the right hand side because if ^ is a flag 
of flats, the negative of a vector w in the weight class of J- can be obtained 
as 

-a; = {-uj\fi vfo) © (-^^Ifz-Fi vfi) © • • • © (-^If^-f,..! vf^-i)- 

To see that the right hand side is contained in —B{M), since vf G —B{M) 
for any flat F, it suffices to show that —B{M) is tropically convex. 

A consequence of the previous description of B{M) is the following. A 
vector X G M" is in —B{M) if and only if, for any r G M, the set (x; r) = 
{i G [n] : j;j > r} is a flat of M. Now take 2;, y G —B{M) and a, 6 G M, and 
let 2: = (a j;) © (6 y), so = min(a + Xi,b + yi). Then, for any r G M, we 
have that (z;r) = (x;r — a) fl (y; r — b). This is a flat in M, because both 
(x; r — a) and (y; r — b) are flats. Thus z G — jB(M). 

Finally, we prove the claim about the vertices of —B{M). If F = Fi n 
• • • n Ffc is an intersection of larger flats, then 

t'F = (0 WFi) © • • • © (0 fFfe) 

SO fi? is not a vertex. Conversely, suppose that we have an equation 

■WF = (ai ■WFi) © • • • © (flfc "WfJ- 

We can assume that Oj 7^ 00 for all i. For each / G F, {vF)f = 00. Thus 
for all i we have {vFi)f = 00; i.e., / G Fj. For each f ^ F, {vF)f = 0. 
Thus for some i we have {vF^f 7^ 00; i.e., / ^ Fj. We conclude that 
F = Fin---nFfc. □ 
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For any tropical polytope P in TP"^^, Develin and Sturmfels j7j gave an 
explicit construction of a nearest point map vrp : TP"~^ — > P, which maps 
every point in tropical projective space to a point in P which is closest to it 
in the ^oo-metric: 

||a;-y||oo= max \{xi - xj) - {yi - yj)\. 

l<i j'<n 

This nearest point map is given as follows. Let x € TP"^^. For each 
vertex v of P, we need to compute A^: the minimum A such that {XQv)®x = 
X. Then 

^p{x) = (XvQv). 

V vertex of P 

Proposition 3. The tropical projection 7r_gj.^^^ maps each vector to to its 
subdominant M -ultrametric: 

^-g{A./)(-'^) = 

Proof. In our case, it is easy to see that = maxj^jj Uf. Therefore 

TT(—Lo)e = min max — w,- 

= —max min ojt, 

eGC* /GC* ^ 

remembering that the cocircuits of M are precisely the complements of its 
hyperplanes. This last expression is the result of applying the blue rule to 
element e. □ 



5 Phylogenetic trees 

Theorem^ is of particular interest when applied to M{Kn), the graphical 
matroid of the complete graph Kn- As shown in the Bergman fan B{Kn) 
can be regarded as a space of phylogenetic trees. Our results thus provide a 
new point of view on a known algorithm in phylogenetics. In fact, it is this 
context that provided the original motivation for our results. 

Let us now review the connection between the Bergman fan B{Kn), 
ultrametrics, and phylogenetic trees. For more information, see |T|. 

Definition. A dissimilarity map is a map 5 : [n] x [n] — > R such that 
S{i,i) = for all i G [n], and S{i,j) = for all i,j € [n]. An ultra- 

metric is a dissimilarity map such that, for all i,j, k G [n], two of the values 
6{i,j),6{j,k) and 5{i,k) are equal and not less than the third. 
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We can think of a dissimilarity map (5 : [n] x [n] ^ R as a weight function 
ios € Rv 2 j on the edges of the complete graph Kn- The connection with 
our study is given by the following result. 

Theorem. ^ A dissimilarity map 5 is an ultrametric if and only if u)s is 
a AI(Kn) -ultrametric. 

As mentioned earlier, the previous theorem is our reason for giving vec- 
tors in B{M) the name of M-ultrametrics. We will slightly abuse notation 
and write 6 instead of lvs; this should cause no confusion. 

Ultrametrics are also in correspondence with a certain kind of phyloge- 
netic tree. Let T be a rooted metric n-tree; that is, a tree with n leaves 
labelled 1, 2, . . . , n, together with a length assigned to each one of its edges. 
For each pair of leaves u, v of the tree, we define the distance dxiu, v) to be 
the length of the unique path joining leaves u and v in T. This gives us a 
distance function dT ■ [n] x [n] R. We are interested in equidistant n-trees. 
These are the rooted metric n-trees such that the leaves are equidistant from 
the root, and the lengths of the interior edges are positive. (For technical 
reasons, the edges incident to a leaf are allowed to have negative lengths.) 




A BCD B 02 c 

Figure 4: An equidistant tree and its distance function. 

Figure 0] shows an example of an equidistant 4-tree, where the distance 
from each leaf to the root is equal to 2. It also shows the corresponding 
distance function, recorded on the edges of the graph K^. This distance 
function is precisely the M(i^4)-ultrametric of Figure ^ This is not a coin- 
cidence, as the following theorem shows. 

Theorem, jlll Theorem 7.2.5] A dissimilarity map J : [n] x [n] is an 

ultrametric if and only if it is the distance function of an equidistant n-tree. 

We can think of equidistant trees as a model for the evolutionary rela- 
tionships between a certain set of species. The various species, represented 
by the leaves, descend from a single root. The descent from the root to a 
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leaf tells us the history of how a particular species branched off from the 
others, until the present day. For more information on the applications of 
this and other similar models, see for example [1] and 

One important problem in phylogenetics is the following: suppose we 
have a way of estimating, in the present day, the pairwise distances ^^(i, j) 
between two species. The goal is to recover the most likely tree T. It is 
well-known how to recover a tree T from its corresponding ultrametric dx 
jlll Theorem 7.2.8]. Of course, in real life, the measured data S{i,j) will not 
be exact. No tree will match the measured distances exactly, and we need 
to find the tree which approximates them most accurately. When proximity 
is measured in the i^o metric: 

\\S -drWoo =max|(5(i,j) -dT{i,j)\, 

Chepoi and Fichet gave a very nice answer to this question, which we now 
review. 

Given a weight function uj on the edges of K^, let 

0Ju{x,y) = min max oj{e). 

paths P from x to y edges e of P 

It is not difficult to see that uju is an ultrametric, known as the suhdominant 
ultrametric of uj. Notice that ojjj = a;^^^"): the formula above is the result 
of applying the red rule to edge xy. 

Write 2e = 1 1 id — W(/||oo — min lu.'(c) — ^(/(e)!, and define a second ultra- 
metric by u}^''{e) = u>u{e) + e for each edge e of Kn- 

Theorem. ]^ Given a dissimilarity map CO on [n\, an t^-optimal ultramet- 
ric for UJ is the ultrametric ujjj"" . 

Our discussion of matroid ultrametrics and tropical projection provides 
a conceptual explanation of the previous theorem. Chepoi and Fichet's goal 
is to construct an -^QQ-closest point to uj in the Bergman fan ISi^Kn). As 
mentioned in Section ^ we can, in fact, essentially solve this problem for 
any tropical polytope via tropical projection. 

In tropical projective space, an £oo-optimal ultrametric is vrgj.^ ^{^), the 

tropical projection of uj onto B{Kn). By Proposition |5J this is also the 
subdominant M(Er„)-ultrametric, uj^^^^"\ Therefore, in real space, an l^- 
optimal ultrametric will be uj^^^"-\ up to a shift by a constant. It is easy 
to see that the optimal constant is — uJ^^^^"^^\oo = £• 
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